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Abstract. We use a si mple a nalytic model to deproject 2-d luminosity functions (LF) of galaxies in the Coma cluster measured 
by Beijersbergen et al. 2002 We demonstrate that the shapes of the LFs change after deprojection. It is therefore essential to 
correct LFs for projection effects. The deprojected LFs of the central area have best-fitting Schechter parameters of M{, = 



-18.31] 
-1.271J 



and au 



-1.27 



+0.018 
0.018' 



-19.79 



+0.14 
0.15 



and ob 



-1.44: 



and M, 



-21.77: 



and a T 



; . The corrections are not significant enough to change the previously observed trend of increasing faint end slopes 
with increasing distance to the cluster center. The weighted U, B, and r band slopes of the deprojected LFs show a slightly 
weaker steepening with increasing projected cluster radius. 



Key words. Galaxies: clusters: individual: Coma (A1656) 
formation. 



Galaxies: luminosity function - Galaxies: evolution - Galaxies: 



1. Introduction 

In the past large field studies of the photometric properties 
of galaxies in the Coma cluster have been based on photo- 
graphic plates (e.g. Godwin et al. 1983; Lugger 1989[ ). More 
recent studies have been made with CCD detectors which pro- 
vide higher photometric accuracy. However, a limitation of 
the use of CCDs has been the relatively small field of view. 
Consequently, these studies have been limited to very sm all ar- 
eas, ranging from ~ 50 arcmin 2 (e.g. Bernstein et al. 1995 



Mo bashe r & Trentham |1998| ) to 1500 arcmin 2 (e.g. Lobo et 
al. 1997), mainly focused on the cluster center. These small 



fields can not be used to study the dependence of the lumi- 
nosity function (LF) on the projected distance from the cluster 
center. Furthermore, lacking the photometric data at larger dis- 
tances from the core it is not possible to correct a central LF for 
projection effects. Accordingly, the published LFs of galaxies 
in the Coma cluster are all projected and are, therefore, lumi- 
nosity distributions rather than actual LFs. It is important to test 
how well these correspond to the actual deprojected LFs. 



Recently, Valotto et al. (2001) used numerical simulations 
to analyse systematic effects in the determination of the galaxy 
luminosity function in clusters. They find that projection ef- 
fects conspire to pr oduce artificially steep faint end slopes. 
Beijersbergen et al. ( 2002| , hereafter BHDH) presented statis- 
tical U, B and r band projected LFs for galaxies in the Coma 
cluster in five annuli out to a projected radius of 1 .4 degrees for 
the B and r bands and 0.74 degrees for the U band. In total an 
area of 5.2 deg 2 in B and r and 1.3 deg 2 in U has been sur- 



veyed. This data offers the interesting possibility to deproject 
actually observed LFs and to obtain true 3-d LFs. 

In this Research Note we use the information of all annuli 
to deproject the 2-d LFs. Our study extends the work in BHDH 
by inferring for the first time deprojected LFs as a function of 
projected distance from the cluster center. In addition, we study 
how the shapes of the LFs are affected by projection by mea- 
suring the faint end slopes as a function of distance to the clus- 
ter center and compare this with results from numerical simu- 
lations. Throughout this Research Note the projected LFs are 
taken from BHDH and are called 2-d LFs. Deprojected LFs are 
named 3-d LFs. 

2. Deprojection of luminosity functions 

Extensive literature exists on the uncovering of three- 
dimensional structure from observations of projected quantities 
(see e.g. Zaroubi et al. 1998 and references therein). Usually 
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methods have been described which were designed for a spe- 
cific type of observation, in particular X-ray measurements. 
Here we use a very simple, analytic geometric deprojection 
scheme (Fabian et al. 198 1 ). This scheme has mainly been used 
for the analysis of X-ray observations, but in principle it is 
suitable for any type of observations. We use this projection 
scheme to deproject the observed 2-d luminosity distributions 
of the galaxies in the Coma cluster to real, 3-d, luminosity func- 
tions. This scheme consists of assuming that the galaxies are 
symmetrically distributed in a sphere. This sphere is divided 
up in a number of concentric shells, in which the galaxy den- 
sity per magnitude is assumed to be constant. The value of the 
3-d luminosity function in a particular shell can then be found 
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The 2-d errors are estimates of the uncertainty in the subtrac- 
tion of the number of field galaxies. These estimates are based 
on very large control fields and therefore very accurate (see 
BHDH for an extensive discussion). The 3-d errors have been 
calculated using Eq. [I]. The corresponding errors are therefore 
dependent on each other and fully determined by the uncer- 
tainty in the subtraction of the number of field galaxies. The 
vertical scales represent volume densities and the 2-d LFs are 
scaled as to match the bright ends of the 3-d LFs. This sim- 
plifies the comparison of their shapes. In annulus V the galaxy 
density is much lower than in the center, but still higher than in 
the field. That is to say, the derived 3-d B and r band LFs still 
suffer from (minor) projection effects due to cluster galaxies 
at radii even larger than 1.7/i^, Mpc. The contamination by 
projection for the 3-d U band LFs is larger since there we miss 
photometric information for radii larger than 0.5/i^qq Mpc. 
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Fig. 1. Outline of the model used to deproject 2-d luminosity 
functions. The figure can be seen as a cut through a sphere and 
as its projection as seen by the observer. The five circles make 



up the projected annuli as used in Beijersbergen et al. 2002. 
Radii are: I: 0-0.2, II: 0.2-0.3, III: 0.3-0.42, IV: 0.42-0.74 and 
V: 0.74-1.4 degrees (1° = 1.22 h^ Q Mpc at Coma distance). 



by subtracting the contributions from the other shells due to 
projection. In Fig. p] we show an outline of this model. 

The number of projected galaxies per magnitude to, Nj )m , 
in annulus i, and the galaxy density per magnitude, 4> m , are 
related as: 



Ni t m = (f>i >m F(R il Ri-i, Ri) + §3,™ 

j=i+i 



[F(Rj, Ri-\, Ri) — F(Rj-i, Ri-i, Ri)] 



(1) 



where R; 
defined by: 



are the radii of the annuli and the function F is 
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With the available 2-d information we can solve for all 
4>i m , i.e. the true 3-d LFs. The N^ m are dependent on the areas 
of the annuli which are accurately known. However, the outer- 
most annuli are not entirely filled by data: annulus V is filled 
for ~ 80% in B and r and annulus IV is filled for ~ 45% in 
U. We omitted annulus IV for the U band and corrected N^, m 
for the B and r bands by assuming a constant projected density 
throughout this annulus. The function F is solely dependent on 
geometry and has values of order 10~ 3 to 10 _1 . Figs. ^| and |] 
show both the 2-d (open symbols) and 3-d (solid symbols) LFs 
for all bands and annuli (all magnitudes are corrected for atmo- 
spheric extinction). The error bars of the 2-d LFs are not shown 
for clarity, but are somewhat smaller than the 3-d error bars. 
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Fig. 4. Power law slopes a as a function of cluster radius. The 
U and r band slopes are offset in radius by -0.03 and +0.03, 
respectively. The faint ends of the LFs become steeper to- 
wards the cluster outskirts. The overall trend is indicated by 
the dashed line. 



3. Discussion 

In the literature LFs of galaxies in clusters are usually given 
in units of number of galaxies per magnitude per area. What 
actually is measured is a luminosity distribution, denoted here 
as 2-d LF. A luminosity function of a sample of galaxies has 
units of galaxies per magnitude per volume. For observations of 
the entire projected volume of a cluster the two have the same 
shape since for all luminosities the same volume is sampled. 
BHDH studied changes in the LF as a function of distance to 
the cluster center. However, they did not correct for the fact 
that fore- and background galaxies, residing in the outer parts 
of the cluster, are projected onto the central cluster regions. 
The question is how much their results have been influenced 



M. Beijersbergen, W. E. Schaap and J. M. van der Hulst: Deprojection of luminosity functions 



3 



so 
o 



-1 



1 1 1 i 1 1 ; i 1 1 1 i '. 
area I 




t — i — i — | — i — i — i — | — i — i — i — | — r 

area I 




J I " " " I I I I I I I L 



t — i — i — | — i — i — i — | — i — i — i — | — r 

: area I 

o9 ■ 



o 




J I I I I I I L 



20 -16 -12 -20 -16 -12 -20 -16 -12 
M r - 5 log h 100 M B - 5 log h 100 M„ - 5 log h 100 

Fig. 2. Two and three dimensional luminosity functions of area I compared. Open symbols represent 2-d luminosity functions 
and solid symbols represent 3-d luminosity functions. The solid lines represent best-fitting Schechter functions to the 3-d data. 
For clarity the error bars of the 2-d luminosity functions are not shown. 
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Fig. 3. Two and three dimensional luminosity functions of areas II to V compared. Open symbols represent 2-d luminosity 
functions and solid symbols represent 3-d luminosity functions. For clarity the error bars of the 2-d luminosity functions are not 
shown. 



by this effect. Below we study how the shapes of the luminosity In the outermost annulus the projected area is converted to 
distributions given in BHDH are affected by projection effects, volume through a factor depending on the adopted geometry 

(radii), i.e. for this annulus the shapes of the 2-d and 3-d LFs 



4 



M. Beijersbergen, W. E. Schaap and J. M. van der Hulst: Deprojection of luminosity functions 



Table 1. Comparison of best-fitting Schechter parameters for 
2-d and 3-d LFs. M* is not corrected for extinction. 



from the 2-d results, it is different at la only and the trend is 
still clearly there. 
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remain identical. From all LFs of the inner annuli the contribu- 
tions to the galaxy counts from surrounding shells is subtracted 
using the model described in Sect. 2. The impact on the result- 
ing LFs can be seen in Figs. ^ and || When the 2-d and 3-d 
LFs are compared binwise the differences, generally, are small. 
Yet, the error bars become larger and it should be kept in mind 
that the scale is logarithmic. Our result can be explained by 
the small values of the F integrals in Eq. |l] combined with the 
density profile of the cluster. In all expressions we find combi- 
nations of relatively high values of the integrals times low den- 
sities, and vice versa, resulting in a minor net correction term. 
However, this correction term is different for each magnitude 
bin causing the shapes of corrected and uncorrected LFs to be 
different. Below, we quantify the differences by studying the 
change in the shapes of the 3-d LFs as a function of position in 
the cluster. The results are compared to the exact same analysis 
for the 2-d LFs. 



3. 1. Dependence of luminosity functions on radial 
distance from the cluster center 

BHDH found that the 2-d central LFs could be well represented 
by Schechter functions. In Fig. || the best-fitting Schechter 
( ]1976| ) j unctions of the 3-d LFs are shown as solid lines. In 
Table [j] we compare the best-fitting Schechter parameters, a 
and M*, for the 2-d and 3-d case. For the r and B bands the 
Schechter fits are again reasonable representations of the data, 
though the best-fitting parameters have changed. The U band 
fit remains a poor representation of the data. 

One of the main results of BHDH was the finding of a sig- 
nificant steepening of the faint end slopes with increasing clus- 
ter radius. A straight line fit {sr + b), where r is the cluster- 
centric radius, to the weighted U, B and r band slopes at each 
radius gives s = 0.11 ± 0.02 Mpc -1 for the 2-d LFs. In or- 
der to compare this with the situation where all 2-d LFs have 
been corrected for projection we fitted power laws (b 10 aM ) 
to the faint ends of the 3-d LFs in Figs. p| and p| . The infalling 
group around NGC 4839 (Neumann et al. 2001) does not have 
a significant influence on any of the derived LFs. The LFs were 
fitted for M r > -20, M B > -20 and Mtj > -19. We omitted 
the slope of the B band LF for area IV since there the LF is not 
well constrained. The results are shown in Fig. Q A straight 
line fit to the weighted slopes of the 3-d LFs at each radius 
gives s — 0.07 ± 0.04 Mpc -1 and is indicated by the dashed 
line. Although this is a somewhat weaker trend than derived 



Valotto et al. (2001 ) used a deep mock catalogue to analyse 
systematic effects in the determination of the galaxy luminos- 
ity function in clusters. They find that projection effects pro- 
duce artificially steep faint end slopes. The cluster core suffers 
most from objects that are projected onto it. The LFs of area I 
shows similar (r and B band) or shallower (U band) faint end 
slopes when they are corrected for projection. We conclude that 
the impact of projection effects on the shapes of the 2-d LFs in 
Coma is not as severe as their simulations suggest. However, 
the 2-d LFs derived for relatively small areas suffer from pro- 
jection effects and should therefore be properly corrected be- 
fore one can reliably interpret their shapes. 

4. Summary and Conclusions 

In this Research Note we have used a simple analytic model 
to deproject 2-d LFs of cluster galaxies. We have applied this 
to the LFs of the galaxies in the Coma cluster measured by 
BHDH. Under the assumptions of a spherically symmetric 
galaxy distribution and a constant galaxy density in each shell 
the 2-d LFs have been corrected for projection effects. The 
main results are summarized by Figs. | and B. Projection ef- 
fects affect the shapes of the 2-d LFs and in general make the 
faint end slopes artificially steeper, although the differences are 
small. The alterations are not significant enough as to arrive at 
different conclusions compared to those inferred from 2-d LFs. 
Furthermore, the previously derived steepening of the faint end 
slopes with increasing cluster radius is only slightly weakened 
when projection effects are taken into account. Overall, we 
have demonstrated by using a simple model that before inter- 
preting shapes of 2-d LFs one needs to correct for projection. 
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